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Abstract 



We study radiative decay of gravitino dark matter with trilinear R-parity viola- 
tions. We show that the branching ratio of the decay of gravitino into monochromatic 
photon can be large enough to explain the observed gamma-ray line from the Galac- 

cn ' tic centre in the Fermi-LAT data without producing too much continuum gamma-ray 

and anti-proton flux. This scenario is realized when the mass of sfermions and the 
trilinear R-parity violating coupling are 0(1 — 10) TeV and 0(10^^ — 10^^) respec- 

K> ■ tively. 



1 Introduction 

Recent studies on the four-year Fermi data have found excess of 130 GeV gamma-ray hne 
from the Galactic Center (GC) [H El El IH El E] . There have been many papers studying 
various possible explanations (instrumental effects [3 El (H [10], pulsar wind effects [11] 
etc.) of this signal but it is most interesting if it is to be interpreted as a dark matter 
(DM) signature [12l [13]. If this interpretation is correct, the 130 GeV gamma-ray line is 
the long-awaited signature of non-gravitational interactions of particle DM. 

In order to explain the 130 GeV gamma-ray line with particle DM, a rather large 
branching ratio of DM annihilating or decaying into monochromatic photon is required, i.e. 
Br(DM — )■ 7) ^ 0.01 p ^ [T5 | [T6]. Otherwise, DM to fermion and gauge boson annihilation 
or decay channels would produce too much continuum gamma-ray and conceal the line 
signal. Moreover, anti-proton flux produced by these channels are constrained by cosmic- 
ray observations jTT] . However, in many cases, the branching ratio of DM annihilating or 
decaying into photons is suppressed because DM does not couple directly to photon. It 
has been shown that for monochromatic photon production channel with standard model 
(SM) particles running in the loop, annihilating DM is typically in tension with the 130 
GeV gamma-ray line scenario |[T8]. 

We consider decaying DM in this letter. Specifically, we study gravitino DM in R- 



parity violating (RPV) supersymmetric (SUSY) models P^ [2U] . v^H With bilinear RPV 
operators, it is difficult to realize the gamma-ray line scenario [H]. We complement the 
previous study by considering trilinear RPV operators. With sfermion masses of O(TeV), 
the tree-level decay rate of gravitino is suppressed and the radiative decay can explain 
the 130 GeV gamma-ray line. Furthermore, there are no overproductions of continuum 
gamma-ray and anti-proton flux. 

The model considered here is consistent with cosmology. The lightest SUSY particle 
of the MSSM (MSSM-LSP) decays into gravitino or other SM particles due to RPV inter- 
actions before the big-bang nucleosynthesis (BBN) begins. This prevents the late decay of 
the MSSM-LSP from spoiling the success of the BBN. We also note that the requirement 
of relatively heavy sfermions does not contradict with the recent discovery of the 126 GeV 
Higgs-like boson at the LHC [221 123]. In fact, negative results on SUSY searches at the 
LHC and a rather heavy Higgs boson favor SUSY models with large sfermion masses. 

The rest of the letter is organized as follows. First, we will discuss the theoretical 
framework of our model. Then, we discuss general aspects of its phenomenology. Next, 
we elaborate how it can explain the 130 GeV gamma-ray line. Before we comment on our 
model and make conclusion, we study cosmological aspects of the model. 



*^ Another viable decaying SUSY DM candidate, axino, has been studied in Ref. [21] . 



2 Gravitino dark matter with R-parity violation 

2.1 Framework 

Let us write down the relevant interaction lagrangian of gravitino: 
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0*,X\ i^u^ A''"^"are sfermion, the corresponding fermion, gravitino and gaugino, respec- 
tively. Dfj^ is the covariant derivative and Fyp " is the field strength tensor. Pl {Pr} is the 
projection operator projecting onto left-handed (right-handed) spinors. All interactions 
are suppressed by the reduced Planck mass Mpi ~ 2.4 x 10^^ GeV . 

Next, we write down the superpotential related to RPV t24j. In the most general form, 
it is 

W = XijkLiLjEk + X'ijj^LiQjDk + X'ljjJJiDjDk + fiiLiHu, (2.2) 

where summation among the indices i,j,k = 1,2,3 which denote the lepton and quark 
generation is implicitly assumed. Lj, Ei, Qi, Di, Ui and H^ are chiral superfields of lepton 
doublet, lepton singlet, quark doublet, down- type quark singlet, up- type quark singlet 
and up- type Higgs doublet, respectively. The first three terms lead to trilinear RPV while 
bilinear RPV arises due to the last term. Xijk, A^^^, Af^^are dimensionless parameters and 
/ij is a parameter with mass dimension one. Hereafter, we will work in the basis where 
^iLiHu is rotated away from the superpotential. This is done by redefining Li and the 
down-type Higgs superfield Hd as L[ = Li — eiHd and H'^ = Hd + eiLi with e^ = /ij//U, 
where fi is the higgsino mass parameter in the MSSM superpotential ^HuHd- Due to 
this redefinition, SUSY-breaking soft terms, including those corresponding to the bilinear 
RPV, 

- Csoft = {BHuHd + BiUHu + m%H,L,Hl + h.c.) + m^Jff^p + m|jL,|2 + ... , (2.3) 

where Li is the scalar component of the chiral superfield Lj, undergo transformation as 
well. In the following, primes of the redefined fields and soft terms are omitted to tidy up 
our notations. In this basis, sneutrinos' vacuum expectation values (VEVs) are typically 
non-zero. By minimizing the scalar potential of sneutrino (including the SUSY-breaking 
soft terms), the VEVs are found to be 



m? TT cos B + Bi sin B 
{h) = - "'^^ ^, (2.4) 

'"'Ui 

tan/3 = Vu/vd is the ratio of the VEVs of the up- and down-type Higgs fields, m?. is the 



sneutrino mass and v = \/v\ -|- v\ ~ 174 GeV. 



V.A;,,<5xlO-'(i::;|i|f) (2.6) 



Note that trilinear RPV terms, LLE and LQD, are also generated by the field redef- 
inition. They are absorbed into the parameters Ajj^and A^ ;,. In summary, we work with 
the superpotential 

W = XijkUL.Ek + KjkLiQjDk + KjkU^D.JDk (2.5) 

and non-zero sneutrino VEVs, which we parametrize as Kj = {vi) /v. 

2.2 Some phenomenological implications and constraints 

In this subsection, we consider general phenomenological aspects of our framework. First, 
we focus on trilinear RPV couplings. Allowing both lepton and baryon number violation 
would lead to proton decay with a very short lifetime. Moreover, as will be discussed in 
detail in the next subsection, gravitino DM's leptonic decays are preferred over hadronic 
decays in order to explain the 130 GeV gamma-ray line. Therefore, we assume baryon 
number conservation by choosing A^'-^ = . 

Another bound on trilinear RPV couplings arises from cosmological considerations. 
Xijk and \[-f^ has to be small enough to prevent wash-out of the baryon asymmetry before 
the electroweak transition. The couplings are, generically j 

Ajg 
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where Msusy is the masses of squarks or sleptons. Other bounds on trilinear RPV cou- 
plings are known to be less stringent [2^ . 

We can no longer distinguish between the lepton doublet and the up-type Higgs dou- 
blet under bilinear RPV. Sneutrino VEVs that mix leptons with gauginos are induced. 
Specifically, neutrinos mix with neutralinos whereas leptons mix with charginos. Bilinear 
RPV's constraints can be deduced from neutrino masses generated by sneutrino VEVs. 
Masses of neutrino are iriy ~ g^{y)'^ jm^^ where m^ is the mass of bino [221 EDI EI]- For 
gaugino masses of 0(TeV),fi;j = (y-^jv < 10~^ is required by experimental bounds on 
neutrino masses. 

In the following, we will study the scenario where trilinear RPV is dominant and 
bilinear RPV is negligible. We note that even if bilinear RPV is absent at tree level at a 
certain energy scale, renormalization group evolution will generate the bilinear terms at 
some other energy scale. Unfortunately, we have not found any model that can naturally 
explain the smallness of bilinear RPV in the literature. On the fiip side, if the 130 GeV 
gamma-ray line can really be interpreted as a signature of gravitino DM with trilinear 
RPV, it should inspire model building efforts towards a theory with such characteristic in 
the future. 

2.3 Trilinear R-parity violation-dominant scenario 

Gravitino undergoes three-body decay via trilinear couplings at the tree level. As will 
be explained in the next paragraph, we will consistently be working with the LLE RPV 



operators. For leptonic decay with an intermediate mass rrirj^^ the decay rate is [20] 

The full analytical result has been worked out in [32]. Gravitino undergoes one- loop 
radiative decay as well |33]. The decay rate scales as 
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when the sferniion masses are large compared to gravitino and lepton masses. The radiative 
decay rate is proportional to the fermion mass. Therefore, RPV couplings involving the 
third generation of fermion gives the largest contribution. We also note that the radiative 
decay is approximately independent of the mass of sfermion running in the loop of the 
decay amplitude. It can be understood from the amplitude of the interaction that involves 
the first two terms of Eq. (12. ip . These terms carry a derivative of the sfermion field that 
brings a momentum flow proportional to the largest loop-mass to the vertex. It balances 



out the contribution of the sfermion mass from the sfermion propagator. |^ Hence, by 
increasing the mass of the sfermion and taking 7/13/2 — 260GeV , one can get the branching 
ratio r('?/^3/2 — ;■ Tt'e) that is large enough to explain the Fermi line. The dependence of 
this branching ratio on the sfermion mass is shown in Fig. [1] 

The apparent 130 GeV gamma-ray line can be explained by a decaying DM of lifetime 
TDMthat satisfies 

rDM/Br(DM ^ 7//) = (1 - 3) x lO^^ sec (2.9) 

and branching ratio Br(DM — )■ 7Z/) > 0.01 pUtj. In this region of parameter, astrophysical 
constraints from diffuse gamma-ray [M] and neutrino spectra [361 ES] ^^^ cilso satisfied. 
As an illustration of our model, we choose the sfermion masses to be m^ ~ 3 TeV. The 
branching ratio of the radiative decay is Br(DM — )■ 7Z/) ~ 0.18. In order to conciliate with 
Eq. (12. 9p . the RPV coupling are needed to be A = A133 ~ 6 x 10"'^. 

We note that a large range of parameters (DM lifetime and branching ratio) is excluded 
by the PAMELA anti-proton data for gravitino DM with bilinear RPV [14j. We do 
not have such concern for gravitino with LLE RPV operators because no anti-proton is 
produced from its decay. 

We now discuss several cosmological constraints on our model. Under trilinear RPV, 
the MSSM-LSP, which we assume to be the bino, may decompose into SM particles via 
tree-level decay. The decay rate is 

5A^a 



"^^See Ref. [33] for a detailed description of the radiative decay amplitude. 
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Figure 1: Dependence of branching ratios of gravitino decay via A133 on the intermediate 
sfermion mass m^. Here, the intermediate sfermions are assumed to be the left-handed and 
right-handed staus. Their masses are degenerate and take the value m^. The decay modes 
are tree-level decay (dashed, red) and radiative decay (blue). The mass of the gravitino is 
m3/2 ~ 260GeV. 



a and Qw are the fine structure constant and the Weinberg angle respectively. We have 
assumed that the masses of left-handed and right-handed stau are degenerate, i.e. m^ = 



"TR 



mv^.The function '?/'(?/) is defined as 

.1/2 
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(2.11) 



For ms/mjj < 10 and m^ ~ ITeV, A has to be greater than 10~^ so that F"^^ < Isec. 
We see that in the region of parameter of interest, BBN is unaffected as bino decays much 
earlier than 1 sec. 

MSSM-LSP decays into gravitino as well. The decay rate is [37] 



F-i 



5xl0^secf^^)'V "^^/^ y 
VlTeV/ V260GeV/ 



(2.12) 



Bino that decays into gravitino contribute to the relic abundance of gravitino. However, 
since F^^ <^ F^^ , only an insignificantly small fraction of bino decays into grav- 
itino. Hence, bino's contribution to the gravitino relic abundance is negligible. 
The thermal relic abundance of gravitino is [5S1 EHl SOI SI] 
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(2.13) 



where Tr is the reheating temperature and rrig is the gluino mass. Since thermal lepto- 
genesis requires Tr > 10^ [42j, our scenario is consistent with thermal leptogenesis for 
m-g ~ O(TeV). 

3 Discussion and Conclusion 

Several comments are in order before we conclude. Gravitino also decays into gauge 
bosons via LLE RPV couplings but branching ratios of these channels are generally equal 
or smaller than that of the radiative decay [53]. This does not alter our results as we still 
obtain branching ratio large enough to explain the Fermi gamma-ray line at large sfermion 
masses without violating astrophysical constraints. 

The morphology of the observed gamma-ray line excess favors annihilating DM but 
decaying DM is acceptable as well [El HI] . More data is required before one confirms or 
rules out the possibility of explaining the gamma-ray line with decaying DM. 

We now briefly discuss the prospect of detection in collider. The lifetime of the MSSM- 
LSP is around 10~^ — 10~^ sec and thus if produced, it will decay outside of the detector. 
Events with missing energy can be recognized as the collider signature. 

In conclusion, we have shown that gravitino dark matter with trilinear R-parity viola- 
tion is capable of explaining the 130 GeV gamma-ray line. Models with the LLE R-parity 
violating coupling are especially advantageous since there is no overproduction of anti- 
proton flux, in contrast with the bilinear R-parity violating scenario. Other astrophysical 
constraints are also satisfied. Furthermore, our model is consistent with cosmology (big- 
bang nucleosynthesis and thermal leptogenesis). The requirement of sfermion masses of 
O(TeV) is well- motivated by the 126 GeV Higgs boson and negative searches for super- 
symmetric particles. 
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